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Abstract. The configuration of the regular magnetic field in M 31 is deduced from radio polarization observations 
at the wavelengths AA 6, 11 and 20 cm. By fitting the observed azimuthal distribution of polarization angles, we 
find that the regular magnetic field, averaged over scales 1-3 kpc, is almost perfectly axisymmetric in the radial 
range 8 to 14 kpc, and follows a spiral pattern with pitch angles of p ~ —19° to p ~ —8°. In the ring between 
6 and 8 kpc a perturbation of the dominant axisymmetric mode may be present, having the azimuthal wave 
number m = 2. A systematic analysis of the observed depolarization allows us to identify the main mechanism 
for wavelength dependent depolarization - Faraday rotation measure gradients arising in a magneto-ionic screen 
above the synchrotron disk. Modelling of the depolarization leads to constraints on the relative scale heights of 
the thermal and synchrotron emitting layers in M31; the thermal layer is found to be up to three times thicker 
than the synchrotron disk. The regular magnetic field must be coherent over a vertical scale at least similar to the 
scale height of the thermal layer, estimated to be h t h — 1 kpc. Faraday effects offer a powerful method to detect 
thick magneto-ionic disks or halosaround spiral galaxies. 
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1. Introduction 

The Andromeda nebula, M31, is the nearest spiral galaxy 
to the Milky Way. Despite its high inclination to the line 
of sight, the large angular size of the galaxy allows de- 
tailed studies of its magnetic field and interstellar medium 
(ISM). In particular, the large scale morphology of the 
magnetic field can be investigated with unmatched pre- 
cision. M31 is thus of prime importance in bringing to- 
gether observational data and theory about galactic mag- 
netic fields. 

Early radio wavelength observations of M 31 at A73 cm 
fPoolev Il969f) and All cm (Berkhuijsen & Wielebinski 
119741 Berkhuiisen ll977|l show the continuum emission con- 
centrated in a ring, at a radius of r ~ 50' ~ 10 kpc. The 
first radio polarization observations at All cm, using the 
100m Effelsberg telescope (Beck et al. 119780 . indicated 
that the magnetic field in the southern part of M31 is 
aligned with the optical spiral arms. Beck 1(1982(1 inter- 
preted the All cm data by comparing the observed po- 
larization angles with a model of the polarized emission 
to reveal a predominantly azimuthal large-scale magnetic 
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field, concentrated in the r ~ 10 kpc 'ring', directed in 
the same direction as the rotation of the galaxy. Faraday 
rotation measures (RMs) from polarization observations 
of the southwestern arm of M31 at AA6, 20 cm confirmed 
the presence of a basically axisymmetric spiral magnetic 
field (Beck et al. 11989(1 . A bisymmetric component of the 
magnetic field was suggested by Sofue & Beck 1(1987(1 from 
an analysis of the deviation of the polarization angles at 
All cm from those expected due to a purely axisymmet- 
ric regular magnetic field; however it is not clear whether 
the inferred bisymmetric mode is statistically significant. 
Ruzmaikin et al. (1990) modelled the All cm polarization 
angles of M31 with an azimuthal Fourier expansion for 
the regular magnetic field and ascertained that deviations 
of the magnetic field from axial symmetry are evident sta- 
tistically and may indicate bisymmetric or higher modes. 
More recently, RMs of 21 background radio sources in the 
field of M 31 were found to be compatible with the same 
magnetic field structure, but extending far away from the 
r ~ 10 kpc 'ring', probably to 5 < r < 25 kpc (Han et al. 
1998). This remains to be substantiated with a statisti- 
cally significant number of sources. 

Recently, Berkhuijsen et al. {ZD03:) presented a new 
A6 cm survey of M 31 and concluded that: the regular com- 
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ponent of the magnetic field is probably as strong as the 
turbulent field; the regular magnetic field has an average 
pitch angle of ~ — f5° in the range 8 < r < 12kpc, with 
a negative value indicating a trailing spiral; gradients in 
Faraday rotation measure may be an important cause of 
depolarization. 

In this paper we seek to take the next logical step in 
understanding the magnetic structure of M 31 by develop- 
ing a detailed and self-consistent description of the mag- 
netic field. We use all of the radio polarization surveys 
(AA6, 11,20 cm) and fit together information on polariza- 
tion angles, Faraday rotation, non-thermal radio emission 
intensities, depolarization and the scale heights of ISM 
components. Our analysis has two main components: de- 
ducing the large-scale geometry of the magnetic field and 
deriving parameters of the magneto-ionic ISM from anal- 
ysis of depolarization of the synchrotron emission. Our 
approach is the latest in a sequence of methods used 
to interpret radio polarization observations of external 
galaxies. Ruzmaikin et al. l|1990fl considered the varia- 
tion of polarization angles at a single wavelength, Sokoloff 
et al. l|1992p extended this approach to multiple wave- 
lengths and Berkhuijsen et al. I|1997|l introduced variation 
in the intrinsic angle of polarized emission in a galaxy. 
We develop a new model, by combining an analysis of 
multi-wavelength polarization angles - based on the ear- 
lier methods - with modelling of the wavelength depen- 
dent depolarization. 

A short description of the data we use is presented in 
Sect. |21 The properties of the synchrotron disk are dis- 
cussed in Sect. |3 In Sect. we use polarization angles 
at AA6, 11,20 cm to deduce the three-dimensional struc- 
ture of the regular magnetic field in M31. The method, 
developed from that used by Berkhuijsen et al. (19511) to 
determine the regular magnetic field of M 51, takes into ac- 
count the intrinsic angle of polarized emission in the disk 
of M31, Faraday rotation by the magneto-ionic medium 
in M 31 and Faraday rotation in the Milky Way. In Sect.[(j] 
we analyze the radial and azimuthal variation in the de- 
polarization between wavelengths A6 cm and A20 cm, and 
derive constraints on the scale heights of the thermal and 
synchrotron emitting disks of M31. This demonstrates a 
new and potentially powerful method for extracting such 
information from radio polarization observations of spiral 
galaxies. A short discussion of the preliminary results was 
presented in Fletcher et al. (2000). 



2. Observational data 

2.1. Radio continuum emission at XX6, 11 and 20cm 

For the analysis in this paper we adopt the follow- 
ing parameters of M31: a distance of 690 kpc (1' = 
196 pc on the major axis), a centre position of RA50 = 
h 40 m l?8, DEC50 = 40°59'46", an inclination angle of 
i = 78° (Braun 1991) where 0° is face-on, and a position 
angle of the northern major axis of 37°. 



Berkhuijsen et al. P3U5|> observed a field of 150' x 70' at 
A6.2cm with the 100m Effelsberg telescope . The original 
resolution was 2.'4. Figure shows the polarized intensity 
smoothed to a beamwidth of 3', along with ellipses show- 
ing the radial range considered in this paper. Preliminary 
results were also discussed by Han et al. (1998) and Beck 

mm . 

The All. lcm map of M31 was obtained with the 
Effelsberg telescope and was published by Beck (1982). 
The original resolution of 4.'4 was smoothed to 5' for this 
paper. The VLA map at A20.5 cm by Beck et al. !|1998t has 
an original resolution of 45", used in the analysis of polar- 
ization angles in Sect.El For the comparison of polarized 
intensities at AA6,20cm presented in Sect. El the A20cm 
map was smoothed to 3', the same as the A6 cm map. Since 
we are interested in wavelength dependent effects (Faraday 
depolarization) we require the same degree of wavelength 
independent depolarization at A6cm and A20cm and so 
the resolutions must be the same (wavelength indepen- 
dent depolarization arises from unresolved fluctuations of 
the polarized emission - see Appendix ^) . 

From the A20cm map in total intensity, smoothed to 
a resolution of 3', and the A6cm total intensity map at 
the same resolution Berkhuijsen et al. (2003) computed a 
spectral index map and maps of thermal and non-thermal 
emission at these wavelengths. Combination of the polar- 
ized and non-thermal emission at each wavelength then 
yields the non-thermal degrees of polarization that are 
analyzed in Sect. |H1 

Missing spacings affect the diffuse emission at A20 cm 
detectable by the VLA. This was corrected in total inten- 
sity with the help of Effelsberg data at the same wave- 
length. A further complication is that M 31 lies behind a 
spur of the Milky Way seen in A20 cm non-thermal emis- 
sion (see., e.g., Grave et al. lTMn Beck et aL HTMll . The 
total emission from the foreground at A20 cm was removed 
using an Effelsberg map of the extended region of sky 
in the direction of M31, but the polarized emission can- 
not be separated into Milky Way and M31 components 
so that a correction for missing spacings was not possi- 
ble in the maps of Stokes Q and U (Beck et al. I1998f) . 
However, strong spatial variation in the Faraday rotation 
intrinsic to M31, shown in Figure [7] for RM derived using 
AA6,ll,20cm data means that at A20cm Stokes Q and 
U originating from M 31 will change rapidly with position 
and hence the effect of missing spacings is probably small 
for the emission from M31. Note that a similar pattern 
is present when RM is determined using only AA6, 11cm 
(Fig. 12 in Berkhuijsen et al. 127)1131) . 

2.2. Data averaging in rings and sectors 

The maps in the Stokes parameters /, Q and U, at each 
of the three wavelengths, were averaged in sectors of 20° 
azimuthal and 2 kpc radial width, in the range 6 < r < 14. 
The size of the sectors was chosen to match the resolution 
of the data at A6 cm. Next we describe how the average 
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M31 6.3cm Effelsberg Polarized Intensity + B-Vectora (%) HPBW=3' 




Fig. 1. Polarized intensity (contours) of M31 with the orientation of the emission B-vector also shown (dashes, not 
corrected for Faraday rotation) with their lengths proportional to the degree of polarization, observed at A6 cm with 
the Effelsberg radio telescope (Berkhuijsen et al. 1201) 33. Note that the foreground RM of — 90radm -2 (Table 2) 
corresponds to Faraday rotation of about 20° so that the intrinsic B-vectors are roughly azimuthal. The beam width 
is 3' and the rms noise is 0.2 mJy/beam. Contour levels are 1, 2, 3, 4, 6 x (5 x 10 -4 ) Jy/beam. A length of B-vectors 
of 3' corresponds to a degree of polarization of 36%. The northern major axis is to the left and the ellipses show the 
radial range of the data analyzed in this paper, 6 < r < 14 kpc. 



Q and U intensities in each sector were combined to give 
the average polarization angle and the average polarized 
emission intensity in each sector. 

2.2.1. Polarization angles 

The polarization angle in a individual sector was calcu- 
lated as ip = i arctan (U)/(Q), where (...) denotes the 
average value of the parameter over the pixels within a 
sector. The resolutions used were 3', 5' and 45" at the 
wavelengths AA6, 11, 20 cm respectively. The errors in po- 
larization angle were computed as the standard deviations, 
within one sector, between all pixels whose intensity is 
stronger than three times the rms noise level. If the num- 
ber of pixels in a sector was below five, the error was 
calculated by averaging several adjacent sectors (this pro- 
cedure was suggested by Berkhuijsen et al. 11997}) . For two 
measurements (both at A20 cm, in the ring 6-8 kpc at 
= 120° and in the ring 8-10 kpc at 9 = 60°) the error 
thus obtained was less than the noise in the maps and 
here the noise error was taken. These average polarization 
angles are analysed in Sect. El 

The A20cm polarized emission from M31 is mixed 
with a substantial amount of emission from the Milky 
Way foreground. At 45" resolution the polarized emission 
from the M 31 'ring' and nucleus is clearly visible and the 
polarization angles are clustered in coherent cells, some- 
times connected with the position of OB associations in 



M31 (see Figs. 2 and 6 of Beck et al. 1998). Thus, the av- 
erage A20 cm polarization angles in sectors with a surface 
area several tens of times larger than the 45" resolution, 
are a reliable measure of the emission from M 31 at this 
wavelength. The foreground Milky Way emission merely 
contributes to the dispersion of angles in a given sector 
and hence to the standard deviation used as our error es- 
timate. 

A further check is applied, by repeating the modelling 
described in Sect. using only the AA6, 11cm data. The 
character of the deduced regular magnetic field does not 
substantially change if the A20 cm is excluded, though nat- 
urally the parameters are less well defined. 

2.2.2. Polarized intensities 

We define the average polarized emission of a sector as 

PI = ((Q) 2 + (U) 2 - l.2(j 2 Q U ) l/2 , where <jq,u is the rms 
noise in Q and U and provides an approximate correction 
for positive bias in PI (Wardle & Kronberg ll974|l . The Q 
and U intensities of all pixels in a sector were averaged to 
compute PI. Errors in non-thermal and polarized intensi- 
ties were estimated as the standard deviation between all 
pixels in a sector as described in Sect. 12.2. Tl 

In Sect. we compare the degree of polarization at 
A20cm, where Faraday effects are strong, with that at 
A6cm, where minimal Faraday rotation occurs. It is nec- 
essary to smooth the A20cm map to the 3' resolution of 
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the A6 cm for this analysis. When smoothed to a reso- 
lution of 3' the 'ring' like polarized emission from M31 
becomes less distinct than at 45" and narrow strips of 
zero polarized intensity become apparent in the A20 cm 
map. These 'canals' are interpreted as depolarization ef- 
fects in the foreground polarized emission of the Milky 
Way by Shukurov & Berkhuijsen ( 2003). The 'contamina- 
tion' of the 3' resolution polarized intensity by Milky Way 
emission is therefore probably more serious than for the 
polarization angles. The azimuthal pattern of the degree 
of polarization at All cm is somewhat similar to that at 
A20cm. However, the difference in the degrees of polar- 
ization at All cm and A6cm is not large enough to allow 
detailed modelling as a check on our results in Sect. H3 
Therefore, the observed A20 cm polarized intensities are 
an upper limit on the emission from M31. 

3. The non-thermal disk 

Our analysis of depolarization in Sect. H3 requires an esti- 
mate of the scale height of the non-thermal disk and the 
discussion of the regular magnetic field, revealed by our 
model in Sect. El is aided by an estimate of the magnetic 
field strengths based on equipartition arguments. In this 
section we derive both of these quantities. 

3.1. The scale height of the non-thermal emission 

In a study of an arm region in the southwest quadrant of 
M 31 Berkhuijsen et al. I|1993jl found that the half-width 
of the arms at A20 cm in the plane of the sky is equal to 
that of the total neutral gas (H1+2H2) suggesting similar 
scale heights for radio continuum emission at A20 cm and 
neutral gas. 

We cannot yet check this for other regions in M31, 
but we can compare the scale heights at A20cm derived 
by Moss et al. I|1998l) with the scale heights of Hi given 
by Braun l|1991[) . Moss et al. (|1998[l determined the scale 
height of the continuum emission from four cuts parallel 
to the minor axis going through the bright 'ring' at about 
20' on either side of the centre. The arms were cut at ra- 
dial distances between 6 and 11 kpc, and the mean of the 
exponential scale heights is 325±43pc. Braun i 1 99 fj de- 
scribed the exponential scale height of the H 1 emission as 
hm = (182 ± 37) + (16 ± 3) r, where the radius r is in kpc 
and /ihi in pc. For the same positions as the radio con- 
tinuum cuts, the mean scale height of Hi is 310 ± 45 pc. 
Hence, the scale height of the radio continuum emission 
at A20cm is the same as that of Hi within errors. At 
this wavelength the width of the radio continuum emis- 
sion is determined by the synchrotron emission, because 
the thermal emission is weak and has a narrower distri- 
bution (Berkhuijsen et al. 2()()()i. Therefore we take the 
synchrotron scale height at A20 cm equal to the H I scale 
height as given by Braun l|1991[) . 

The synchrotron scale height depends on frequency as 
iy-°- 25 , as observed in NGC891 (Hummel et al.ESB and 
M31 (Berkhuijsen et al. I1991J) . thus the scale heights at 



All cm and A6cm are somewhat smaller than at A20cm 
(see Table HJ). 

3.2. The equipartition magnetic field strength 

The transverse component of the total field strength -Btot_i_ 
(the quadratic sum of the regular and turbulent compo- 
nents) can be evaluated from the intensity of the non- 
thermal emission assuming, for example, equipartition be- 
tween the energy densities of magnetic field and cosmic 
rays (see Pacholczyk 1970, Longair I1994fl . However, we 
use a fixed integration interval in cosmic-ray energy rather 
than a fixed interval in radio frequency (see Beck et al. 
1996). In this case, and for a non-thermal spectral index 
a n ~ 1, the equipartition field strength is identical to the 
minimum-energy field strength. The polarized intensity 
yields the strength of the transverse regular field, B±; the 
transverse turbulent field strength b± is then found from 
b\ = B% ot± — B\ . The values of B tot and the regular field 
strength B one obtains by deprojection assuming that B 

1 /2 

is oriented parallel to the plane of M 31, and b = (| bj) 
assuming statistical isotropy. As Faraday effects are small 
at A6 cm, we evaluated the field strengths from the A6 cm 
data. 

In Table we show the average equipartition field 
strengths in four 2 kpc- wide rings covering the bright emis- 
sion from M31 between 6 kpc and 14 kpc radius. We used 
a non-thermal spectral index a n = 1 (Berkhuijsen et al. 
2003 ) , and the standard ratio of relativistic proton to elec- 
tron energy density k — 100. The line of sight through the 
emission layer was taken as L = 2he/ cosz with i = 78°; 
we note again that the synchrotron scale height depends 
on A. 

The magnetic field strengths derived only weakly de- 
pend on the errors in /, L and k (as the power l/(3+a n ) ~ 
1/4). The main uncertainties are in L and k (about 50% 
each) so that the uncertainty in the derived field strengths 
in Table [1] is about 20%. 

In each ring we also calculated the average magnetic 
field strength in the sectors described in Sect. 12.21 within 
the errors -Btot, B and b are constant in azimuth. 



4. Overview of the method 

A short overview of the method we use may help the reader 
follow the main part of the paper. We develop two linked 
models in the following two sections. First, an analysis of 
the average polarization angles is used to deduce the un- 
derlying structure of the regular magnetic field in M31. 
One of the parameters in this model, £a in Eq. Q, can be 
estimated from a second model of the Faraday depolar- 
ization. However the second model, of the depolarization, 
uses rotation measures derived in the first model. We will 
try to find solutions that satisfy both models and are con- 
sistent with each other i.e. the parameter £\ is the same, 
within errors, in each model. 
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-Dtot 


13 


/ 




he 




(kpc) 


(mjy/beam) 


(%) 


(AiG) 


(pG) 


(jmG) 


(pc) 


(PC) 


6- 8 


3.72±0.06 


30±1 


7.3 


4.9 


5.4 


220 


290 


8-10 


4.71±0.05 


33±1 


7.5 


5.2 


5.4 


240 


330 


10-12 


4.19±0.05 


31±1 


7.1 


4.9 


5.2 


270 


360 


12-14 


2.71±0.05 


35±1 


6.3 


4.6 


4.3 


290 


390 



Notes: r is the radial range, Ie the average non-thermal radio intensity per beam area at A6cm (HPBW=3'), Pe the non- 
thermal degree of polarization at A6cm obtained from the average polarized intensity per beam area divided by Ie ; he and /120 
the exponential scale heights of the synchrotron emission at A6cm and A20cm, respectively; Btot, B and b are the average 
equipartition strengths of the total, regular and turbulent fields, respectively. The uncertainty in the derived field strengths is 
about 20%. See Sect.|3for further explanations. 



5. The 3D structure of the regular magnetic field 

In this section, we deduce the regular magnetic field in 
M 31 from polarization angles of synchrotron emission at 
AA6, 11,20 cm. The method used is an extension of that 
employed by Berkhuijsen et al. Q1997[) and is only briefly 
described here. 



5.1. The model 

The polarization angle of synchrotron emission is given by 

ijj = ijj ( B) + A 2 RMi( B) + A 2 RMfg, (1) 

where tpo is the intrinsic polarization angle, RMi is the 
Faraday rotation measure in the galaxy, RMf g is the 
Faraday rotation in the Milky Way and A is the wave- 
length. 

The cylindrical components of B = (B r , Bg, B z ) are 
expanded in Fourier series in the azimuthal angle 9, 

B r = B sinpo + Bi sinpi cos(6* — 

+ B 2 sinp 2 cos2(6> - /3 2 ), 
Bg = B cospo + Bi cospi cos(6> — f}\) (2) 

+ B 2 cos p 2 cos 2(6* - /3 2 ), 
B z = B z0 + B zl cos{0 - p zl ) + B z2 cos2{6 - (3 z2 ), 

where B m and B zm are the amplitude of the mode with 
azimuthal wave number m in the horizontal and vertical 
fields, p m is the pitch angle of the m'th horizontal Fourier 
mode (i.e. the angle between the field and the local cir- 
cumference) and /3 m and f3 zm are the azimuths where the 
non-axisymmetric modes are maximum. Expressions for 
^0 and RM m terms of the expansions shown in Eq. (T2J 
are given in Eqs. (A3) and (6) of Berkhuijsen et al. ( 1997). 
We note here that ipo depends on magnetic field compo- 
nents in the sky plane whereas RM depends on those along 
the line of sight. Therefore, fitting Eq. allows us to 
obtain all three components of B. Since the observed po- 
larization angle depends on RM, i.e. on the product of 
the magnetic field strength, thermal electron density and 
the path length, the amplitudes of the Fourier modes are 



obtained from fitting in terms of the variables R ri 
dimension is radrn" 2 : 



0.81 



B„ 



( n e) 

1 cm _: 



L 

1 pc 



whose 



(3) 



where (n e ) and L are the average density of thermal elec- 
trons and the line of sight path length through the thermal 
disk in a given ring. 

Only a fraction of the synchrotron emitting disk may 
be visible at a given wavelength due to Faraday depo- 
larization. Therefore, observations of polarized emission 
at different wavelengths probe the galactic disk to dif- 
ferent depths and our analysis can reveal variations in 
the disk parameters along the line of sight. This has al- 
lowed Berkhuijsen et al. I1997[) to reveal a two-component 
magneto- ionic structure in M51 comprising a disk and a 
halo. M31 does not have an extensive synchrotron halo 
(Grave et al. I1981f) . and so we consider one-component 
(i.e. disk only) fits where the galactic disk is probed to 
different depths at different wavelengths. Correspondingly, 
Faraday rotation is scaled by a wavelength dependent fac- 
tor < 1, so 



RM; = £xRM d , 



(4) 



where RM^ is the Faraday rotation measure produced 
through the whole disk thickness (observable at short 
wavelengths), and t;\ can be understood as the fraction 
of the disk thickness transparent to polarized emission at 
a wavelength A. In Sect. El we discuss depolarizing mech- 
anisms in detail, and from models of the observed depo- 
larization we adopt ^6 = 1-0, £11 = 0.96 and £20 = 0.75 
(see Eqs. IS] and El. The fitted parameters describing the 
magnetic field were found to be rather insensitive to the 
adopted values of £>, . 

Using Eqs. (JTJ, © and Q we fit the modelled, three- 
dimensional B to the observed polarization angles in a 
ring, simultaneously for all wavelengths, by minimizing 
the residual 



(5) 
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Table 2. Parameters of the fitted model and their 2a 
errors. RMf g is the Faraday rotation measure arising in 
the Milky Way, R m and p m are the amplitude and pitch 
angle of the mode with wave number m, and p) m is the 
azimuth where a mode with azimuthal wave number m 
is maximum. The minimum value of the residual and the 
value of x 2 are shown for each fit in the bottom lines. 



Units Radial range (kpc) 







6-8 


8-10 


10-12 


12-14 


RMf g 


rad m 2 


-93 ±5 


-99 ±5 


-93 ±5 


-89 ±4 


Ro 


rad m -2 


+83 ±7 


+96 ±9 


+ 115 ±9 


+99 ±6 


Pa 


deg 


-13 ±4 


-19±3 


-11 ±3 


-8 ±3 


R2 


rad m -2 


+45 ±10 








P2 


deg 


-2 ±12 










deg 


-43 ±7 








s 




58 


59 


62 


62 


x 2 




63 


63 


65 


65 



where ip n is the observed angle of polarization, ip{9 n ) 
the modelled angle and a n are the observational errors. 
The x 2 test is used to ensure that the fit, for all wave- 
lengths in a ring, is sufficiently close to the measured 
angles. Application of the Fisher test verified that the 
fits are equally good at each individual wavelength (see 
Bcrkhuijsen et al. H997f) . This model is aimed at analysis of 
a global structure of the magnetic field and is not devised 
to capture local details in the field structure. Therefore, 
where a few data points deviate very strongly from the 
general pattern they can be discarded to obtain a statis- 
tically good fit. The number of points discarded for this 
reason was 6% of the total measurements and most of 
these points occur where tp varies very strongly with 9, 
leading to underestimated errors. The exclusion of points 
is discussed further in Sect. 15.31 

We determine the errors in the fit parameters by vary- 
ing them independently and in paired combinations to 
determine the parameter ranges consistent with the \ 2 
test. For fits requiring a small number of parameters, we 
checked these error estimates by plotting contours of the 
residual S in the parameter space. The resulting errors, 
quoted below, are all 2a deviations. 

5.2. Results of fitting 

Figures [3 to [S] show the variation of observed polarization 
angles (tp, measured anti-clockwise from the local radial 
direction in the plane of M31) with azimuthal angle 9 
and the fits for each ring. The fitted parameters are given 
in Table [5] Generally we find that an axisymmetric field, 
lying parallel to the galactic midplane provides the best 
fit to the measured polarization angles. For the innermost 
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100 200 300 

Azimuthal angle, 8 (°) 



Fig. 2. Polarization angles (ip, measured from the local 
radial direction in the plane of M31) against azimuth (9) 
for the ring 6-8 kpc. Fit (solid) and observations (squares 
with error bars, horizontal lines with error bars show ex- 
cluded points) are shown for A6cm, All cm and A20.5cm 
from top to bottom. The error bars show the la devia- 
tions. 

ring a weaker, 7r-periodic (m = 2) mode is added to the 
dominant axisymmetric (m = 0) mode. The m = 2 mode 
will produce a 7r/2 periodicity in RM. 

The fitted RMf g is constant, within errors, between 
adjacent rings and varies weakly across the whole radial 
range in agreement with the expected small fluctuations 
in foreground RM from our Galaxy in the direction of 
M31 (Han et al. 1998|). This is an important reliability 
check for the model; the values of RMf g in Table [21 were 
independently derived for each ring by fitting a non-linear 
model to the observational data. It is reassuring that there 
is agreement between rings within errors and with earlier 
estimates. The value of RMf g is broadly consistent with 
earlier estimates of -88 ±2 rad m" 2 (Beck Egg) , -100 ± 
33radm~ 2 (Ruzmaikin et al. HMty . -93±3radm" 2 (Han 
et al. 11353)1 and -92 ±3 rad m" 2 (Berkhuijsen et al- EDDSjl - 

The median amplitude of the axisymmetric mode i?o 
reaches a maximum at R ~ 11 kpc, the radius of the well 
known bright radio 'ring' of M31. However, the maximum 
is only marginally pronounced, and the values of Rq only 
show radial variation at the 2a level. This implies that 
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Fig. 3. As in Fig. but for the ring 8-10 kpc. 

the synchrotron ring in M 31 is prominent either because 
the synchrotron emissivity depends on a high power of the 
magnetic field strength or because the density of relativis- 
tic electrons is higher in the ring. The underlying max- 
imum in the magnetic field itself is very weak, of about 
20%, or even less if the thermal electron density has a 
maximum in the ring. 

We now describe the fits for each ring in detail. 

5.2.1. The ring r = 6-8 kpc 

A combination of a strong m — mode perturbed by a 
weaker m — 2 mode provides a good fit for the innermost 
ring (i.e., one that satisfies both the x 2 an d Fisher tests). 
The rotation measure in this ring varies by a factor of 3 
between the maxima (~ 130radm~ 2 at = 130°, 310°) 
and minima (~ 45radnT 2 at 6 = 50°, 230°). If (n c )L 
is about constant in the ring, the field strength varies by 
the same factor. The pitch angle of the m = 2 mode is 
small but leads to a variation of ±10° in the mean pitch 
angle of the regular magnetic field, p — arctan(B r / Be), 
with minimum pitch angles of p — —25° and maxima of 
p = -9° at 6 = 40, 220° and 140, 320°, respectively. 

To achieve this fit we excluded two data points (at 
A6cm the sectors 8 — 0° and 160°) out of 54. Both sec- 
tors are in the region of a vary rapid change in tp, so it is 
plausible that the error in ip is underestimated in the two 
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Fig. 4. As in Fig. but for the ring 10-12 kpc. 



sectors. If we try to obtain a good fit for the combination 
m = + 1, it is necessary to exclude four measurements 
(at A6cm 6 = 0°, 160°, 340° and at All cm 9 = 180°) 
and the fitted RMf g = — 113±2radm -2 is not consistent 
with RMfg in the other rings. Nine measurements must 
be excluded in order to achieve a fit using only the ax- 
isymmetric m = mode, so the addition of three extra 
parameters describing the m — 2 mode is supported by 
the use of seven extra data points. 

A possible explanation for the m = 2 mode in this 
ring can be that the disk inclination angle i is different 
from that in the other rings. Braun l|1991f) argues that the 
inclination angle of the H I disk varies significantly along 
radius in M 31. 

Another, more plausible possibility is that the m = 2 
component is a response to the two armed spiral pattern, 
but restricted to the thin magneto-ionic disk. The latter 
restriction is needed to explain why this magnetic field 
model does not deliver a good fit to the Faraday depo- 
larization in this ring discussed in Sect. l(T2l As we argue 
there, the depolarization, due to a Faraday screen, occurs 
in the upper layers where the field is basically axisymmet- 
ric 
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Fig. 5. As in Fig.Elbut for the ring 12-14 kpc. 



5.2.2. The rings 8-10, 10-12 and 12-14 kpc 

A satisfactory fit using only the axisymmetric m = 
mode is found for each of these rings. The mode ampli- 
tude reaches a weak maximum in the ring 10 — 12 kpc and 
then decreases in the outermost ring. The pitch angle of 
the regular magnetic field becomes smaller (i.e., the field 
becomes more tightly wound with increasing radius (see 
Sect. Oil. 

The fit for the ring 8—10 kpc requires the omission 
of four measured ip out of the total of 54, two near the 
major axis at 9 = 340° for AA6, 11 cm, and two at A20cm, 
9 = 200°, 220°. For the ring 10-12 kpc three measurements 
must be omitted to achieve the m = fit, two on the major 
axis (at A6cm, 9 = 160° and All cm, 9 = 0°), along with 
the sector 9 = 40° at A20cm. Finally, the measurements 
at A6 cm 9 = 0°, 340° and A20 cm 9 = 340° are omitted in 
the outermost ring. 

Figure shows a face-on view of the galactic disk with 
the sector grid and the fitted regular magnetic field vectors 
shown in each sector. The azimuthal component of the 
field is stronger than the radial component in all sectors 
(that is, the pitch angle is rather small). The effect of 
the 7r-periodic, m = 2, mode in the innermost ring can 
be clearly seen in the varying length and direction of the 
magnetic field vectors. 




Fig. 6. Face-on view of M31 showing sectors and regu- 
lar magnetic field vectors obtained from the fits shown in 
Table The grid radii are 6, 8, 10, 12 and 14 kpc. The 
length of the vectors is proportional to B. 



5.3. Excluded measurements 

In order to include all of the observations in any of the 
rings, we find that more than two extra modes must be 
added to the magnetic field models discussed above. For 
example, in the ring 10-12 kpc we cannot achieve a good 
fit with the combination m = + 1 + 2 even though an 
extra 6 parameters are used to try and accommodate three 
previously excluded measurements. This strongly suggests 
that either (i) the excluded sectors are not dominated by 
any large-scale structure but by localised perturbations of 
a more regular underlying pattern or (ii) the errors in the 
omitted polarization angles are underestimated. 

All but two of the excluded measurements lie close to 
the major axis of M 31. Here, the detected polarized emis- 
sion is weakest as the small pitch angle of the regular mag- 
netic field means that its component perpendicular to the 
line of sight, B±, is small near the major axis. Also, near 
the major axis of a highly inclined galaxy with a strongly 
azimuthal regular magnetic field, the observed polariza- 
tion angle (in the sky plane) changes rapidly. These ef- 
fects can lead to underestimation of the errors in tp for 
sectors near the major axis. Furthermore, any deviation 
from an axisymmetric field (e.g. due to inter-arm bridges) 
near the major axis of M 31 contributes to the line-of-sight 
magnetic field and distorts the smooth pattern of Faraday 
rotation measures. 



5.4. Magnetic pitch angles 

The pitch angles of the regular magnetic field are p ~ —17° 
between 6 < r < 10 kpc and then become smaller with in- 
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creasing radius, reaching p ~ —8° in the ring 12 < r < 
14kpc. These values are more reliable than earlier esti- 
mates - more data are used in the modelling and interpre- 
tation methods have improved - but are in broad agree- 
ment with the results of Beck (1982), Ruzmaikin et al. 
( 1990) and Berkhuijsen et al. 1|2003[1 . The regular magnetic 
fields maintained by galactic dynamo action must have a 
non-zero pitch angle, since the dynamo generates both ra- 
dial and azimuthal magnetic field components (Shukurov 
2000}. The sign, magnitude and radial trend of the mag- 
netic field pitch angles are in broad agreement with the 
predictions of a range of dynamo models for M31 ({2000 ). 

Observations of CO (Guelin et al. l2000|) and Hi (Braun 
1991) have been fitted with logarithmic spirals tracing the 
gaseous arms with a constant pitch angle of ~ —7°. In 
those nearby spiral galaxies where density waves are be- 
lieved to be present, the regular magnetic fields generally 
follow the spiral structure (see Beck l|1996(l and references 
therein). The difference between the magnetic and spiral 
arm pitch angles for 6 < r < 12 kpc may be because den- 
sity waves are absent or very weak in M31. A detailed 
comparison with the spiral structure, seen e.g. in the CO 
line emission, is required to clarify the relation between 
the magnetic and gas spirals. 

6. Depolarization 

The observed degree of polarization of non-thermal emis- 
sion from external galaxies is generally less than the in- 
trinsic maximum of Po — 0.75 for a completely regular 
magnetic field structure. The reduction in the degree of 
polarization can be due to the physical properties of the 
ISM in the galaxy and to effects arising from the finite 
size of the telescope beam. By investigating depolariza- 
tion mechanisms we can recover information about the 
ISM. 

A convenient measure of depolarization is the ratio of 
relative polarized intensities at two wavelength, i.e., 



DP 



Ai/A 2 



(6) 



where P(X) is the degree of polarization at a wave- 
length A; DP^/^ = 1 means no depolarization between 
the two wavelengths. A variety of depolarization mecha- 
nisms in radio sources are discussed by, e.g., Burn 



Pacholczyk (1970) and, specifically for spiral galaxies, by 
Sokoloff et al. (1998). A description of several concur- 
rent depolarization mechanisms can be rather compli- 
cated. Wavelength-independent depolarization and that 
due to Faraday rotation (and so wavelength-dependent) 
can be easily isolated as they result in independent fac- 
tors in the total depolarization; in Eq. © the wavelength 
independent contributions to depolarization at Ai and A2 
are equal to each other and cancel, so that the observed 
DP is a measure of depolarization due to Faraday ef- 
fects. Among wavelength dependent depolarization mech- 
anisms, depolarization in a Faraday screen and in the syn- 
chrotron source can be disentangled because they occur 



in non-overlapping regions. However, distinct Faraday de- 
polarization mechanisms that occur within the same vol- 
ume cannot be represented as independent factors in the 
total depolarization in the general case. An approxima- 
tion that allows one to separate the internal Faraday dis- 
persion from other depolarizing effects in the synchrotron 
source has been suggested by Sokoloff et al. I|19981 their 
Sect. 6.3) and called the 'opaque layer approximation'. In 
Appendix El we briefly describe the different depolariza- 
tion mechanisms affecting observations of external galax- 
ies and give the equations used later in this section to 
model the observed depolarization. 

Depolarization of the non-thermal emission must be 
carefully considered when interpreting the data. For ex- 
ample, the synchrotron disk can be transparent to polar- 
ized emission at short wavelengths but opaque at longer 
wavelengths (see e.g. Berkhuijsen et al. 199||. Therefore, 
the amount of Faraday rotation is no longer proportional 
to A 2 . (This is the motivation for introducing the param- 
eter £>, in Sect. El) First though, we look at the observed 
depolarization in a qualitative way. Then we attempt to 
construct a model for the observations, in terms of param- 
eters describing the state of the ISM. 

6.1. The dominant depolarization mechanism in M 31 

The ring 10 < r < 12 kpc is chosen for an initial, closer 
look at depolarization. In Fig. we show the azimuthal 
variation of some key properties in this ring. RM and B±_ 
have been derived from the polarization angle model pre- 
sented in Sect. B± has been obtained assuming that 
(n c )L (see Eq. 3) is constant in azimuth, and normalized. 
We also show the observed degrees of polarized emission 
at A6cm and A20cm (P 6 and P20 respectively). The pat- 
tern of RM versus azimuthal angle is determined by the 
geometrical variation of Bu , with the strongest RM near 
the major axis where the regular magnetic field lies along 
the line of sight to M31. Note also that the sine-like vari- 
ation of RM results in the strongest gradients in RM 
lying near the minor axis. Furthermore, Berkhuijsen et 
al. ( 2003) noted that the azimuthal variation of the polar- 
ized emission at A6 cm is almost completely due to the geo- 
metrical variation of B± with azimuthal angle. FigurcsQ) 
and [3; clearly show that this also holds for Pq, with the 
Pq highest near the minor axis where B± is strongest. In 
contrast, the degree of polarization at A20cm, P20, has a 
less marked azimuthal variation. If anything, the pattern 
of P20 is the inverse of Pq, but with a lower amplitude. The 
wavelength dependent depolarization, DP 2 o/6j is obtained 
by dividing P 20 (Fig. 01) by P 6 (Fig. Et). 

Can we recognize the signature of any of the depolar- 
ization mechanisms discussed in Appendix ^ in the ob- 
servations? Before considering a model based on a combi- 
nation of effects it is instructive to consider each of these 
mechanisms separately. 

The effect of wavelength independent depolarization is 
removed by considering Eq. JHJ). By comparing the ratio of 
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Fig. 7. For the ring 10 < r < 12kpc. (a) Rotation mea- 
sures derived from the fitted magnetic field shown in 
Tabic [3 (b) The normalized amplitude of the component 
of the regular magnetic field lying in the sky plane, -Bj_, 
derived from the fitted magnetic field shown in Tabic [2] 
(c) The observed degree of polarization (P&— polarized in- 
tensity/total nonthermal intensity) at A6cm. (d) The ob- 
served degree of polarization at A20cm (f^o)- 



the observed degrees of polarization at A20 cm and A6 cm, 
DP20/6 = P20/P6, with that expected from Eqs. I|A.1|) to 
l|A.5|) . we identify which wavelength dependent depolar- 
ization mechanisms are dominant. 

The observed DP 20 / 6 , plotted in Fig. [HJ has a marked 
azimuthal variation with strong depolarization at A20 cm 



near the minor axis 



90° and 270°, where DP 



20/6 



0.1) and less depolarization on the major axis (9 = 0° 
and 180° where DP 20/6 > 0.4). Note that the observed 
DP 2 o/6(#) is roughly proportional to the derivative of 
RM(6>) (Fig. Q . 

Gradients in the foreground RM due to magnetic fields 
in the Milky Way, RMf g , in the direction of M31, are 
weak (Han et al. I1998|l and unlikely to cause the observed 
variation of RM and DP with azimuth in M31. Thus, de- 
polarization must occur within M31. For the rest of the 
analysis of depolarization we consider the RM intrinsic to 
M31, RMi = RM - RM fg . 

The smooth, sinusoidal azimuthal variation of RM 
(Fig. Hi) can be completely accounted for by an azimuthal 




100 200 300 

Azimuthal angle, 6 (°) 

Fig. 8. Observed (squares with error bars) and expected 
(solid line) depolarization between A20 cm and A6 cm, 
DP20/6: f° r the ring 10 < r < 12kpc assuming vari- 
ous depolarization mechanisms. The azimuthal angle 9 is 
measured counterclockwise from the northern major axis. 
Depolarization due to (a) differential Faraday rotation us- 
ing Eq. I|A.1J| , (b) RM gradients in the synchrotron source 
and differential Faraday rotation from Eq. (|A.4|> and (c) 
RM gradients in a Faraday screen given by Eq. lfO)l arc 
represented by solid lines. RM and the gradient in RM are 
derived from the fitted magnetic field described in Sect.|5| 
We conclude that foreground RM gradients, illustrated in 
panel (c), dominate other wavelength dependent depolar- 
ization mechanisms. 



variation of B\\ deduced in Sect.|Sl indicating that {n c )L is 
indeed roughly constant in azimuth. The turbulent mag- 
netic field, 6, derived using the equipartition approach de- 
scribed in Sect. is also constant in azimuth for each ring. 
Therefore the dispersion in RM, ctrm, and hence depolar- 
ization due to Faraday dispersion, is roughly constant at a 
given radius, and the azimuthal variation in DP 2 o/6 can- 
not be explained by Faraday dispersion [Eqs. (|A.2|I and 
(|A.3|l ]. This does not mean that Faraday dispersion is in- 
effective in M31, but rather that the strong azimuthal 
pattern in DP 20 /6 cannot be explained by this mechanism. 

The remaining wavelength dependent depolarizing 
mechanisms are all caused by the regular magnetic field: 
differential Faraday rotation, RM gradients within the 
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emitting layer, and RM gradients in a foreground Faraday 
screen. The first two effects are unavoidable while the third 
effect requires the existence of a 'thick disk' of magnetic 
fields and thermal gas invisible in synchrotron emission. 
Differential Faraday rotation in the source is strongest 
near the major axis where the line-of-sight magnetic field 
-By is maximum, resulting in a depolarization pattern very 
different from that observed (Fig. [Ht) • The azimuthal gra- 
dient in _B|| is also maximum near the minor axis (where it 
changes sign). Therefore, depolarization due to gradients 
in RM in the synchrotron source is strong near the mi- 
nor axis and weak near the major axis, but still does not 
overcome the differential Faraday rotation that produces 
a different pattern (Fig. ISJs) . On the other hand a fore- 
ground Faraday screen does not produce any differential 
Faraday rotation, and so depolarization due to the RM 
gradients in a foreground screen is dominant, producing a 
correct pattern shown in Fig. 

Thus, the global pattern of the azimuthal variation of 
DP 2 o/6 can only be reproduced by depolarization due to 
RM gradients in a Faraday screen (the bottom frame of 
Fig. |HJ. This mechanism must be the dominant cause of 
the azimuthal pattern in wavelength dependent depolar- 
ization. This is true in the whole radial range 6 < r < 
14kpc. Berkhuijsen et al. (2003) found that contours of 
RM and DP n / 6 are often perpendicular to each other 
where they cross (see their Fig. 14) and noted that this 
suggests RM gradients as an important cause of depo- 
larization. Earlier, Berkhuijsen & Beck (19iIQJ) found that 
RM gradients were primarily responsible for depolariza- 
tion in the southwestern quadrant of M31 and Horcllou 
et al. I| 19921 observed that contours of DP and RM are 
perpendicular at crossing points for the galaxy M51. 





1.0 




0.8 


a 

° 


0.6 


Q_ 




O 


0.4 




0.2 




0.0 




1.0 




0.8 


a 

s 


0.6 


Q. 







0.4 




0.2 




0.0 




1.0 




0.8 


■ 


0.6 






O 


0.4 




0.2 




0.0 




1.0 




0.8 




0.6 


0- 




O 


0.4 




0.2 




0.0 



(0) q=0.1 DP,„=0.30 




JL 



S 



100 



200 300 

i i 1 

(b) q=0.4 DP,„=0.30 




_ I I ' ' ' 






1 1 _ 

(c) q=0.7 DP in =0.30 I 






1 




Li , 




f 






h 







100 



200 



300 



(d) q=1.0 DP,„=0.30 




100 200 
Azimuthal angle, 8 (°) 



300 



Fig. 9. Observed (squares with error bars) and expected 
(solid line) depolarization between A20 cm and A6 cm, 
DP20/6: due to RM gradients and differential Faraday ro- 
tation using Eqs. J7J|, (fTT|) and (fT^|l . for the ring 10 < 
r < 12 kpc and with different ratios of the synchrotron to 



The minima in DP 20 /6 produced by the Faraday screen thermal disk scale heights, q = h syn /h 



20/6 

are noticeably deeper than those observed (at 9 ss 90° 
and 270° in Fig. IHfc) ■ As discussed in Sect. 16.21 this can be 
explained by other, less important depolarization mecha- 
nisms. 



6.2. The thermal and synchrotron disk scale heights 

We have identified RM gradients in a Faraday screen as 
the dominant depolarizing mechanism responsible for the 
observed azimuthal pattern of DP 2 o/6 in M31. The fit to 
observations in Fig. [HJ; can be improved by including other 
depolarizing effects, especially Faraday dispersion. Also, 
the effectiveness of the Faraday screen depends upon its 
relative thickness, compared to that of the synchrotron 
emitting layer. Now we attempt to recover information 
about the relative heights of the emitting and Faraday 
rotating layers from fitting the depolarization. 

A full description of depolarization due to the regu- 
lar magnetic field (i.e., differential Faraday rotation, RM 
gradients inside the synchrotron emitting layer and in a 
Faraday screen) is given by the product of Eqs. I|A.4|) and 
l|A.5(l . The intrinsic Faraday rotation measure RMj and 
its increment A RM across each sector can be calculated 



tii - 



from the fits for the regular magnetic field discussed in 
Sect. GD We split RMi into two components, RMd aris- 
ing within the synchrotron disk and RM S arising in the 
part of the thermal layer above the synchrotron disk (see 
Fig. I10|) ; the scale height of the synchrotron layer is taken 
from Sect. 13. ll The first component will produce depolar- 
ization due to differential Faraday rotation, but the latter 
will only contribute to Faraday screen effects. The gradi- 
ent in RMj is similarly split into A RMd and A RM S . In 
terms of these variables, the degree of polarization with 
allowance for Faraday dispersion, differential Faraday ro- 
tation and rotation measure gradients in both the thermal 
disk and Faraday screen is given by 



P(A) = P 0j Pin(A)Pex(A) exp 2iRM 0s A 2 -2 (ARM S A 



1 ,2 x sin(2A RMdA 2 s, , 

exp(4i RM 0d A s) — f n — - ds 







2ARM d A 2 s 



(7) 



where subscript zero refers to a value at the sector cen- 
tre. We can reasonably assume that depolarization due to 
Faraday dispersion in the emitting layer is much stronger 
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Fig. 10. Sketch showing the scale heights of the thermal 
disk, /i t h, synchrotron disk, h syn , and the depth in the 
thermal disk from which polarized emission is visible at 
A20cm, Az. The galactic midplane is at z — 0, i.e., at the 
bottom of the figure. The deeper layers of the synchrotron 
disk are invisible in polarized emission because of inter- 
nal Faraday dispersion. This sketch is only an illustration: 
note that Az can exceed h t h if more than half the disk 
thickness is visible. However, h t h — h syn < Az < h t h + h syn 
with the extreme values corresponding to DP 2 o/6 = and 
1, respectively. 



than Faraday dispersion in the foreground screen (tur- 
bulent magnetic fields and thermal gas density will be 
stronger near the mid-plane) and so Pi n P C x ~ -Pin- 

Since the galactic disk may be opaque to polarized 
emission at longer wavelengths, mainly due to internal 
Faraday dispersion, the effective path length can differ 
from that suggested by the disk scale height (c.f. Sokoloff 
et al. 119981 Berkhuijsen et al. I1997|l . We use the 'opaque 
layer' approximation of Sokoloff et al. (119981 Sect. 6.3) 
to describe the visible depth Az in terms of the depolar- 
ization due to internal Faraday dispersion assuming that 
all the observed polarized emission at A20 cm arises from 
an upper layer in the synchrotron disk. Figure I1UI shows 
how hth, h syn and Az are related. The path lengths over 
which the observed polarized emission is produced are 
Az — (hth — h syn ) at A20cm (here Az is a function of 
A) and 2ft, syn at A6 cm, where the disk is assumed to be 
transparent to polarized emission. Then a crude estimate 
of Az in terms of the observed degrees of polarization fol- 
lows from assuming that depolarization due to internal 
Faraday dispersion is constant for all sectors in a ring: 



A:. 



^th 



h syn (2DP in -l) , 



(8) 



where DP m = P in (20cm)/P in (6cm) is depolarization due 
to internal Faraday dispersion alone (see Sect. 3.3.3 in 
Berkhuijsen et al. 11997(1 . Since some depolarization due to 
other mechanisms occurs within Az, this yields minimum 
values for the thickness of the visible layer. Equation (JHJ 
then gives the parameter £\ in Eq. J3J, in terms of q = 
h s 



t 



"syn/^th and PP; n , via (Berkhuijsen et al. H997[) 

Az — h s 



A 



\ ' 1 



'syn 



Hh 



so we have 

a = 1 + q(DP ir 



1)- 



(9) 



(10) 



Table 3. Parameters derived from the model of Faraday 
depolarization and their 2a errors. 



Ring (kpc) 


q — hgyn/kth 




fhh (kpc) 


6-8 


No statistically good fit found 


8-10 


0.6 ±0.1 


0.1 ±0.1 


0.55 ±0.10 


10-12 


n 4+ - 4 


o.3tS:I 


0.90i° ;*° 


12-14 


o s +0 - 5 


0.2 ± 0.2 


-i q+2.6 
1 -°-0.8 



The parameter £\ links the modelled regular magnetic 
field described in Sect. from which we obtain RM 
and A RM, and the model for Faraday depolarization de- 
scribed in this section. Our aim is to obtain satisfactory 
fits for both models using the same £\ in each. 

For a thermal layer thicker than the synchrotron disk 
(the configuration that produces a foreground Faraday 
screen), Az > h t h — h syn , Berkhuijsen et al. (1997) showed 
that 

RM 0d = \ RM AZ - ( ^ th - fesyn) = q RMo DP in , (11) 



RM 0s = RM 



hth hs y n 



>lh 



(1 - q) RMo, 



(12) 



where the final equalities result from substitution of 
Eq. (JSJ. Using Eqs. l|TTjl and (JT^J we can express P(X) in 
Eq. as a function of q = h syn /hth, the ratio of the scale 
heights of the synchrotron and thermal disks, and DP m . 
We fit the values of q and PPi n by comparing DP 2 o/e ob- 
tained from Eq. (JJJ to the observed values in each sector. 
Figure compares the calculated DP 20 /6, f° r Q between 
0.1 and 0.4, with observed values in the ring 10-12 kpc, 
and shows how increasing the relative scale height of the 
synchrotron disk reduces the effect of the Faraday screen 
and enhances depolarization due to differential Faraday 
rotation, which is strongest on the major axis. 

Figure HP shows the azimuthal variation of the depo- 
larization for the best fitting q for each ring. The fitted 
values of q — h syn /h t h and DP m are given in Tabic 
where the quoted 2a errors represent the extent of the 
(q, PPi n )-parameter space within the relevant x 2 contour. 
These errors are large, but it is remarkable that we can 
successfully model Faraday depolarization in such a com- 
plex system using such a simple, two parameter model. 
Despite the uncertainty in the precise values of the model 
parameters, the key result of this section is robust; the 
strong azimuthal pattern of depolarization can only be ex- 
plained by a Faraday screen acting within M 31 and hence 
the thermal electron layer must be significantly thicker 
than the synchrotron emitting layer. 

For 8 < r < 14 kpc the modelled DP 20 /6 reproduce 
the observations well. For each of these rings the fitted 
DP 2 o/6 meets the \ 2 test for statistical significance at the 
2a level. In the two rings at largest radii the results of 
the depolarization modelling are fully consistent with the 
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Fig. 11. Observed (squares with error bars representing 
ltr error) and expected (solid line) depolarization between 
A20cm and A6cm, DP 2 o/6, due to Faraday effects using 
Eq. J7J|, and 1|12|) . for each of four radial rings. The ex- 
pected depolarization is fitted by varying the ratio of the 
synchrotron to thermal disk scale heights, q = ft-syn/^th 
and depolarization due to internal Faraday dispersion, 
DP m . The fitted values of q are 0.6, 0.6, 0.4, and 0.3 and 
of DP m 0.1, 0.1, 0.3 and 0.2 in panels (a), (b), (c) and 
(d), respectively. The solid line was obtained using RM 
derived from the regular magnetic field fits of Table [2] 



polarization angle model used to deduce the regular mag- 
netic field structure in Sect. 03 The two models are linked 
by the parameter £\ - a weighting for the depth in the 
emission layer visible at long wavelengths - in Eqs. 
and l |10| l. For the rings 10-12 kpc and 12-14 kpc, q and 
DP hl in Table give £>o = 0.7 ± 0.3 and f 2 o = 0.8+ - 2 , 
respectively. In Sect. we adopted £20 = 0.75 for all of 
the rings and discrepancies of order ±0.1 have a negligi- 
ble effect on the fitted magnetic field parameters given in 
Table |2 

For the ring 8-10 kpc, q and DP in give £ 20 = 0.5 ± 0.2, 
whereas the best fit to the polarization angles in Table |2 
requires £20 = 0.75. We can achieve self consistency be- 
tween the magnetic field and depolarization models by dis- 
carding more measured polarization angles in Sect. [5] i.e. 
by making the model of the magnetic field worse. However, 
the main problem with the depolarization model in this 



ring is that around the north end of the major axis our 
method of averaging the data gives zero average polar- 
ized emission at A20 cm in three sectors. (The polarized 
intensity is averaged from maps smoothed to 3' resolu- 
tion; the polarization angles are derived from 45" maps 
at A20 cm and do not suffer from this problem.) Without 
measurements at both ends of the major axis the fitted 
DP20/6 favours a model with stronger Faraday dispersion 
i.e. a model that has a less prominent double minimum. 
For these reasons we prefer to retain the magnetic field 
model obtained with £ = 0.75 - to keep £20 the same for 
each ring - and accept that the depolarization model for 
this ring is poorer than for 10-12 kpc and 12-14 kpc. 

The quality of the fit is clearly bad in the ring 6- 
8 kpc (Fig. ITTk 'l where the magnetic field, deduced in 
Sect. contains both the axisymmetric (m = 0) and 
the quadrisymmetric (to = 2) components as given in 
Table [21 In Sect. 15.2.11 we show that this ring may have 
a more complicated regular magnetic field structure than 
the purely azimuthal fields in the other rings. The mod- 
elled azimuthal patterns of RM and the gradient in RM 
are rather complicated in the ring 6-8 kpc and no good fit 
can be obtained. The results would be better if we used a 
simpler fit involving a purely axisymmetric magnetic field. 
However, as explained in Sect. l5.2~Tl nine polarization an- 
gle measurements must be discarded to make an m = 
magnetic field model, and the consequent degrading of the 
regular magnetic field model is not justified. 

Using h syn , at A20cm, from Tabled we can estimate 
h t h from the fitted values for q = h syn /h t h- The results are 
shown in Table These scale heights are about a factor 
of two or more greater than previously expected in M31, 
where the low star formation rate and absence of a radio 
halo were thought to imply the likely absence of a thick 
ionized disk (Walterbos & Braun lT^Mjl , 

We emphasize that the gradients in rotation measure 
producing most of the depolarization in M31 are due to 
the highly axisymmetric regular magnetic field that we 
find from an analysis of polarization angles in Sect. |SJ 
For simplicity, in modelling the depolarization we assumed 
that the regular magnetic field has the same configuration 
and strength throughout the full vertical extent of the 
thermal layer (including the synchrotron emitting disk). If 
the regular magnetic field strength or the thermal electron 
density has a maximum above the emitting disk (i.e., at 
z > 300 pc), the RM required to produce the observed 
depolarization can be generated in a thinner layer and hth 
will be lower than estimated above, but still /ith > ^syn- 

In Sect. 12.2.21 the limitations of the A20 cm polariza- 
tion data when smoothed to 3' were discussed; foreground 
emission from the Milky Way cannot be subtracted from 
the emission from M31 and so the A20cm polarized in- 
tensities are upper limits. The values of q = h syn /hth and 
DP m shown in Tabled were derived assuming that all of 
the polarized emission at A20cm comes from M31 and so 
are upper limits on q and DP m . 

The corresponding lower limits (i.e., giving depolar- 
ization stronger than required) can be obtained assuming 
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that the emission from M 31 is nearly completely depo- 
larized at A20cm. Without Faraday depolarization, the 
A20 cm polarized emission from M 31 will have the same 
azimuthal pattern as the A6cm PI shown in Fig. |7f c). 
Total depolarization (DP 2 o/6 = 0) will occur when DP lri 
and differential rotation are just strong enough to depolar- 
ize the emission on the major axis (8 = 0°, 180°) and RM 
gradients are just sufficient to depolarize emission from 
the minor axis (9 = 90°, 270°). From Fig. Eta) we estimate 
that, for the ring 10< r <12 kpc, complete depolarization 
will occur if q ~ 0.1 and DP ln ~ 0.1. These are the lower 
limits on q and DP in . 

The regular magnetic field must be coherent in z over 
at least the scale height of the thermal disk, and we have 
shown that the latter must exceed that of the synchrotron 
disc. This poses the intriguing question of why the cosmic 
rays in MSI are confined to a layer several times thinner 
than the regular magnetic field. One possible answer re- 
lies on the usual assumption of equipartition between the 
cosmic ray and magnetic field energy densities. Then the 
synchrotron emissivity depends upon the fourth power of 
the magnetic field and so Kb ~ ^h syn ~ 1.5 kpc. This scale 
height is in good agreement with h t h derived from our 
analysis of depolarization (at least for the two rings with 
the most reliable model of DP 2 o/6)- I n M31, the magnetic 
field is well ordered with B ~ b and there is no signifi- 
cant vertical component of the magnetic field (see Sect.EJ- 
This may be sufficient to suppress diffusion of cosmic rays 
perpendicular to the disk plane and so constrain them to 
the same layer as their sources. 

6.3. Thermal electron densities 

Using Eq. © and the equipartition regular magnetic field 
strengths given in Table ^ the rotation measures from 
Table El and the thermal disk scale heights of Table El 
we can derive average thermal electron densities for M31 
in the radial range 8 < r < 14 kpc. This gives (n c ) ~ 
0.008,0.007 and 0.004 cm- 3 for the rings 8-10, 10-12 and 
12-14 kpc, respectively. These values refer to the upper 
layers of the thermal electron layer, z <; h syn ~ 200- 
300 pc, that act as the Faraday screen. 

Electron density closer to the midplane can be ob- 
tained from the amount of depolarization due to Faraday 
dispersion between A20cm and A6cm, DP ln ~ 0.1 
as obtained above. Using Eq. ljA.21) with b — 5 /iG, 
L = 200 pc and d = 50 pc, we obtain ctrm — 
550 ((n^^/lcm -3 ) radm~ 2 and then DP in = 0.1 cor- 
responds to n c ~ 0.1 cm~ 3 . This estimate is compatible 
with that obtained by Walterbos & Braun jX§94) from H a 
emission measures of the diffuse ionised gas, (n c ) ~ 0.08- 
0.04 cm -3 with a filling factor 0.2. 

Thus, the equipartition magnetic field strength, rota- 
tion measures and the scale heights of the thermal disk 
derived in our models produce an estimate for (n c ) that 
is in broad agreement with (n c ) obtained from completely 
different data and methods. 



Berkhuijsen et al. (2003) note that there is little cor- 
relation between RM and thermal emission in M31 and 
suggest that the small filling factor of Hn regions may 
be the reason. This is consistent with our conclusion that 
much of the Faraday rotation in M31 is produced in a 
Faraday screen. 

6.4. Review of the method 

In fitting the modelled to observed polarization angles in 
Sect. El we use the parameter to account for the par- 
tial opacity of the galaxy's disk to polarized emission at 
A20 cm. In order to estimate £a we need to know the ratio 
of the scale heights of the synchrotron and thermal disks, 
q = h sya / hth, but the values for q deduced in Sect. 16.21 
make use of RM calculated from the fits of Sect. El We 
used an iterative approach to try to obtain a model con- 
sistent with both the observed depolarization and polar- 
ization angles. This method was successful for the two 
outer rings, after one iteration, but not for the rings 6— 
8 kpc and 8-10 kpc. For these rings we adopted £ = 0.75 
from the self-consistent models of the rings 10-12 kpc and 
12-14 kpc. 

7. Summary 

Sensitive, high resolution, multi-wavelength radio polar- 
ization observations have been used to study the magnetic 
field of M 31, between the radii of 6 and 14 kpc. The pow- 
erful method of using polarization angles to uncover the 
regular magnetic field structure was supplemented by a 
systematic analysis of depolarization to produce a model 
of the regular magnetic field which is consistent with all 
of the radio polarization data for 10 < r < 14 kpc. 
Our main conclusions are as follows : 

1. The regular magnetic field in M31 is axisymmetric 
to a very good approximation. 

2. The magnetic field has a significant radial compo- 
nent at all radii and so is definitely not purely azimuthal. 
Vector lines of the regular magnetic field in the radial 
range 6 < r < 10 kpc can be approximated by trailing 
logarithmic spirals, with the pitch angle p ~ —16°. The 
magnetic spiral becomes tighter at large radii, with \p\ de- 
creasing to p ~ —7° at r = 12-14 kpc. The magnitude and 
trend of magnetic pitch angles is in broad agreement with 
those expected from dynamo theory (Shukurov 2000). 

3. Analysis of the azimuthal pattern of the wavelength 
dependent depolarization reveals that a Faraday active 
screen lies above the synchrotron emitting disk of M31. 
The diffuse thermal disk is thicker than previously ex- 
pected, with a scale height of h%h ~ 1 kpc. 

4. The scale height of the regular magnetic field is at 
least equal to hth- 

5. The magnetic field in M31 extends inside and out- 
side of r ~ 10 kpc, as found by Han et al. ( 1998), and does 
not have a strong maximum at this radius. The bright ra- 
dio ring is a result of a high density of cosmic ray electrons. 
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6. The equipartition field strengths are about 5 fiG for 
both the regular and turbulent field components, without 
significant variation between 6 kpc and 14 kpc radius. 

7. Faraday rotation measures and equipartition field 
strengths are in agreement for average electron densities 
of 0.008-0.004 cm -3 . The electron densities inferred from 
the Faraday dispersion measures are ~ 0.1cm -3 , close 
to the average electron densities found by Walterbos & 
Braun Ijl994|) . This suggests that the diffuse ionised gas 
is mainly responsible for the Faraday rotation, with little 
contribution to RM from Hn regions. 

Our analysis of depolarization in M 31 is the most ex- 
tensive undertaken to date for a spiral galaxy, and shows 
that the theory of radio depolarization developed by Burn 
(1966) and Sokoloff et al. (1998 ) can be used not only to 
identify the causes of depolarization, but also to reveal 
properties of the diffuse ISM in external galaxies. 

Acknowledgements. We thank Marita Krause for comments 
following careful reading of the manuscript. 

Our results are based on observations with the Effelsberg 
100m telescope of the MPIfR. 

AF was funded by a PPARC studentship at the University 
of Newcastle, where much of this work was undertaken. 
Financial support from NATO (Grant PST.CLG 974737), 
PPARC (Grant PPA/G/S/2000/00528) and a University of 
Newcastle Small Grant are gratefully acknowledged. 

Appendix A: Depolarization mechanisms 




Fig. A.l. Sketch showing wavelength independent depo- 
larization. The double headed arrows represent E-vectors. 

Wavelength independent depolarization is 

caused by tangling of magnetic field lines in the emitting 
region (Figure IA.1J) . The intrinsic polarization angle 
of synchrotron radiation is perpendicular to the local 
magnetic field orientation and so tangled magnetic field 
lines result in emission at a range of polarization angles 
within a single telescope beam. As long as the beam sizes 
are equal, the degree of depolarization due to tangled 
magnetic field lines will be the same at all wavelengths. 

Faraday rotation by both regular and turbulent mag- 
netic fields results in wavelength dependent depolar- 
ization. It is useful to consider separately Faraday effects 
within the synchrotron emitting layer and Faraday rota- 
tion in regions where there is no emission, i.e. within a 
Faraday screen (Figure 
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Fig. A. 2. Sketch of synchrotron emitting layer and a fore- 
ground Faraday screen. 
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Fig. A. 3. Differential Faraday rotation occurs within the 
synchrotron emitting layer. Emission from different depths 
along the same line of sight undergoes different amounts of 
Faraday rotation, causing depolarization. For clarity the 
sketch separates emission from different depths. 



The regular field in the synchrotron emitting layer 
causes depolarization by differential Faraday rota- 
tion, whereby polarized emission from different depths 
along the line of sight is rotated by different amounts 
(Figure IA.3|) . In a slab with uniform magnetic field and 
electron density the degree of polarization is (Burn 1966 
Sokoloff et al. ITMsT) 



Pr, 



Po 



sin(2RMA 2 ) 



2RMA 2 



(A.l) 



where RM = 0.81(n e ) B\\L/2 is the observed rotation 
measure in units of radm -2 , with (n e ) the average ther- 
mal electron density in cm -3 , B\\ the component of B 
parallel to the line of sight in fiG and L the path length 
through the Faraday active emitting layer, in pc. 

The presence of unresolved, turbulent magnetic field 
means that polarized emission along different lines of sight 
within the telescope beam undergoes different amounts 
of Faraday rotation (Figure IA.4(I . When the emitting 
and rotating layers coincide, the effect is called internal 
Faraday dispersion, and the degree of polarization is 
given by Sokoloff et al. ( 1995J) as 



Pn 



1 



exp(-2a| M A 4 ) 



2cr RM^ 4 



(A.2) 
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Fig. A. 4. Sketch illustrating Faraday dispersion. When 
this effect occurs within the synchrotron emitting layer 
it is called internal Faraday dispersion; occurrence in 
a Faraday screen (no synchrotron emission) is external 
Faraday dispersion. Note that the required variation of 
rotation measure within the beam could also be caused 
by fluctuations of thermal electron density. In that case 
the magnetic field could be totally regular but Faraday 
dispersion would still occur. 



where crj| M = 0.81(n e b) 2 2Ld with d the correlation scale 
(half the turbulent cell size) of the turbulent magnetic 

mm 



field in parsecs. Burn (1966) gives the depolarization due 
to external Faraday dispersion (due to turbulent fields 
in front of the synchrotron source) as 



P cx = P exp(-2cr^ M A 4 ). 



(A.3) 




magnetic field 



Fig. A. 5. Depolarization due to gradients in rotation 
measure: if RM varies across the beam, different parts of 
the beam area contribute with different polarization an- 
gles. Note that variations of thermal electron density or 
magnetic field strength can also produce RM gradients 
even in a totally regular magnetic field. 



Gradients in rotation measure across the beam 
cause depolarization that is especially strong when the 
resolution of observations is low (Figure rO)l . For de- 
polarization by RM gradients within the synchrotron 



source, including the effect of differential Faraday rota- 
tion, Sokoloff ct al. (1998) obtained 



-PAin — Pq 



exp(4z RM A 2 s 



sin(2 ARMA 2 s) 
2 A RM A 2 s 



ds 



,(A.4) 



where RMo is RM at the centre of the beam, ARM is 
the increment in RM across the beam and the normalized 
integration variable s describes the line of sight within 
the synchrotron disk, < s < 1. Equation (|A.4|> holds for 
both resolved and unresolved RM gradients (Sokoloff ct 
al. I1998|l . Depolarization by a resolved RM gradient 
in a Faraday screen is given by Sokoloff et al. ( 1995J) as 



-Pagx — Pq 



exp 



2iRM A 2 - 2 (ARM A 2 ) 5 



(A.5) 
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